o 

=5 



Systematic study of isoscaling behavior in projectile 
fragmentation by the statistical abrasion-ablation 
model 

§ '. D. Q. Fan^l , Y. G. Ma, C. Zhong, C. W. Ma, X. Z. Cai, J. G. 

Chen, W. Guo, Q. M. Su, W. D. Tian, K. Wang, T. Z. Yan, W. 
Q. Shen 

, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 

Cn ' 201800, People's Republic of China 

{Sj . Abstract. The isospin effect and isoscaling behavior in projectile fragmentation have 

been systematically investigated by a modified statistical abrasion-ablation (SAA) 
■ model. The normalized peak differences and reduced isoscaling parameters are found 

, to decrease with (-Zp ro j — Z)/Zp TO j or the excitation energy per nucleon and have no 

significant dependence on the size of reaction systems. Assuming a Fermi-gas behavior, 
the excitation energy dependence of the symmetry energy coefficients are tentatively 
, extracted from a and (3 which looks consistent with the experimental data. It is 

pointed out that the reduced isoscaling parameters can be used as an observable to 
study excitation extent of system and asymmetric nuclear equation of state in heavy 
ion collisions. 
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x. 

1. Introduction 



The process of projectile fragmentation has been studied extensively for investigation 
of reaction mechanisms in heavy ion collisions at intermediate and high energies 
[HOmmOinilT]. It is also one of the most important methods to produce exotic 
nuclei. Recent advances in experiments using radioactive ion beams with large neutron 
or proton excess have lead to the discovery of halo structure |8j [9]. Since then 
interest in the study of very neutron-rich and proton-rich nuclei has grown due to 
their anomalous structures. In addition, the studies of isospin physics have become 
a very popular subject. The isospin effects of various physical phenomena, such as 
multifragmentation, flow, pre-equilibrium nucleon emission, etc., have been extensively 
reported [TU1 [HI 121 H3J [HI HSJ [HI [T7]. The studies have shown that isospin effect 
exists in nuclear reactions induced by exotic nuclei but it may disappear under certain 
conditions. Our previous calculations by using the modified statistical abrasion-ablation 
(SAA) model have demonstrated that the fragment isotopic distribution shifts toward 
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the neutron-rich side for neutron-rich projectile, but the shift decreases with the increase 
of the parameter (Zp ro j-Z)/Zp ro j or the violence of nuclear reaction. This isospin effect 
of fragmentation reaction on the fragment isotopic distribution will disappear when 
(Zp ro j-Z)/Zp ro j becomes larger than 0.5 [T8| [T9]. 

Recently, study of the nuclear symmetry energy has become a very important 
topic in nuclear physics. It is well known that the nuclear symmetry energy is very 
significant for investigation of the nuclear equation of state and a variety of astrophysical 
phenomena. The isoscaling approach for light fragment composition produced in the 
multifragmentation of very hot source has become an important method in heavy 
ion collisions since it can isolate the nuclear symmetry energy in the fragment yields 
[20~1 |2"H [221 [23]. The scaling law relates ratios of isotope yields measured in two different 
nuclear reactions, 1 and 2, B,2i(N,Z) = ^(iV, Z)/Y\(N, Z). In multifragmentation 
events, such ratios are shown to obey an exponential dependence on the neutron number 
iV or proton number Z of the isotopes or isotones characterized by three parameters a, 
P and C [20]: 

R2i(N, Z) = = Cexp(aN + (3Z), (1) 

here C is an overall normalization constant. In the grand-canonical limit, a and (3 will 
have the form, 

a= ^H [( |) 2 -(|)^1^5lA[(|n (2) 

and 

Z- N- 

where Csym is symmetry energy coefficient (MeV), (^-) 2 or (^) 2 (i=l,2) means the 
square of charge or neutron number over mass number for system 1 and 2. T is the 
temperature of the system in MeV. This behavior is attributed to the difference of isospin 
asymmetry between two reaction systems in similar nuclear temperature. Since the 
symmetry energy determines nuclear structure of neutron-rich or neutron-deficient rare 
isotopes, studies on the isoscaling behavior can be used to probe the isospin dependent 
nuclear equation of state [I1IH1E21E31E1E51I1H 

So far, the isoscaling behavior has been studied experimentally and theoretically for 
different reaction mechanisms. However, most studies focus on the isoscaling behaviors 
for light particles with Z=2— 8. A few studies on the heavy projectile-like residues in 
deep elastic collisions and fission fragments have been reported [29], [301 EH E21 [331 Ell- 
in this paper, we will present our studies on systematic behaviors of the isospin 
effect as well as isoscaling features for projectile-like fragments in the framework of 
statistical abrasion-ablation model. Extraction of the symmetry energy coefficient from 
the isoscaling parameters will also be investigated. 



Systematic study of isoscaling behavior in projectile fragmentation 



3 



2. Model description 

The statistical abrasion-ablation model can describe the isotopic distribution well [7j. 
In the SAA model , the nuclear reaction is described as two stages which occur in two 
distinctly different time scales. The first abrasion stage is fragmentation reaction which 
describes the production of the pre-fragment with certain amount excitation energy 
through the independent nucleon-nucleon collisions in the overlap zone of the colliding 
nuclei. The collisions are described by a picture of interacting tubes. Assuming a 
binomial distribution for the absorbed projectile neutrons and protons in the interaction 
of a specific pair of tubes, the distributions of the total abraded neutrons and protons are 
determined. For an infinitesimal tube in the projectile, the transmission probabilities 
for neutrons (protons) at a given impact parameter b are calculated by [7j 

W~ b) = exp{-[Al(r - b)a nk + (r - 6)a pk ]}, (4) 

rp rp rpi 

where D is thickness function of the target, which is normalized by / d rD^ = N 
and / d 2 rDp = with and referring to the neutron and proton number in 

the target respectively, the vectors r and b are defined in the plane perpendicular to 
beam, and oVu is the free nucleon-nucleon cross sections (k', k=n for neutron and k', 
k=p for proton). The thickness function of the target is given by 

rp /" + 00 

d£(t)= dzp k ((r 2 + z 2 y/ 2 ), (5) 

J — oo 

with being the neutron (proton) density distribution of the target. So the average 
abraded mass at a given impact parameter b is calculated by the expression 

(AA(b)) = Jd 2 rDl(r)[l-t n (f-b)] 

+ Jd 2 rDl(r)[l-t v (f-b)]. 

The second stage is the evaporation process in which the system reorganizes due to 
excitation. It deexcites and thermalizes by the cascade evaporation of light particles 
using the conventional statistical model [35]. The excitation energy for projectile 
spectator is estimated by a simple relation of E* = 13. 3( A A(b)) MeV where 13.3 
is a mean excitation energy for an abraded nucleon from the initial projectile. This 
excitation energy was given by the statistical hole-energy model as described in Ref. [35J. 
After the evaporation stage, we can obtain the final fragments which are comparable 
to the experimental data. By introducing in-medium nucleon-nucleon cross section and 
optimizing computational method given in Ref. (36J, [181 EE E7], it can give a good 
agreement with the experimental isotopic distributions [HJ HH [37] . Comparison of the 
SAA model calculations with the experimental isotopic distributions for Z = 30—32 from 
44 AMeV 86 Kr+ 27 Al [38] and Z = 43 - 45 from 790 AMeV 129 Xe+ 27 Al [39] is shown in 
Fig. [H The results shown in this figure and all the following figures are referring the final 
fragments after the evaporation stage. For 86 Kr, all isotopic distributions are normalized 
by the same factor in order to compare with the experimental yields. For 129 Xe, the 
calculated production cross sections are compared with the data directly. From this 
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Figure 1. Comparison of isotopic distributions between the SAA model and the 
data. The isotopic distributions for selected charge numbers from 44 AMeV 86 Kr+ 27 Al 
(upper panel) and 790 AMeV 129 Xe+ 27 Al (lower panel). The dots are the experimental 
data taken from Ref.[38 for 44 AMeV 86 Kr+ 27 Al and Ref.09] for 790 AMeV 
129 Xe+ 27 Al, the lines are the results calculated by the SAA model. 




Figure 2. The normalized peak difference AApg^/ AAp ro j of the fragment isotopic 
distribution as a function of (^p ro j — Z)/Z-p TO y Different symbols are used for 
projectiles with different charge number as shown in the legend. 



figure, we can see that the SAA model can reproduce the experimental data both at 
intermediate and high energies quite well. The isospin effect and its disappearance in 
projectile fragmentation for m ^ m Ar at intermediate energies have been predicted by this 
model and confirmed by the experimental data [T9] . 

3. Calculations and discussion 

The model predicts that strong isospin effect exists in the isotopic distributions produced 
by projectiles with same charge number but different mass number [T8t [T9] . In order 
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Figure 3. Yield ratios R.2i(N, Z) of fragments from the reactions of 124/,112 Sn+ 112 Sn 
at 6QA MeV versus N for the selected isotopes (left panel) and Z for the selected 
isotones (right panel) . Different symbols are used for different isotopes and isotones as 
shown in the legend. The lines represent the exponential fits. For details see text. 



to do a systematic study of the isospin effect in projectile fragmentation, reactions 

Of 4 0/36 Ar> 48/40 Ca; 64/58^ 86/78^ 124/112g n and 129/136 Xe Qn 112g n ftt q QA Me y are 

simulated by the SAA model. Since the isotopic distributions from two projectiles have 
similar shape but a shift in mass, their peak positions will be one of the most sensitive 
quantities for the isospin effect. Thus we extract the peak position by Gaussian fit to 
the fragment isotopic distribution for each charge number Z as in Ref.[T8l [19] . The 
normalized differences of the peak position from two projectiles AAp ea jj./ AAp ro j as a 
function of (Zp ro j — Z) / Z-p TO j are shown in Fig. [2j Here AAp ro j is the mass number 
difference between the two projectiles with same charge number and AAp ea k is the 
peak position difference of the fragment isotopic distribution produced by these two 
projectiles. Zp ro j and Z are the charge number of the projectile and the produced 
isotopes. AAp ea k / AAp ro j exponentially decreases as the increase of (Zp ro j —Z)/ Zp ro j 
which is same as our previous conclusions [19] . The dependence of AAp ea k/AAp ro j on 
(Zp ro j — Z)/Zp ro j shows a very slight difference among different size projectiles. 

To study systematic behaviors of the isoscaling phenomena, the yield ratios 
R,2i(N,Z) are made using the convention that index 2 refers to the more neutron- 
rich system and index 1 to the less neutron-rich one. As an example, Fig. [3] shows 
the yield ratios R 2 i(N,Z) as a function of neutron number N for selected isotopes 
and Z for selected isotones from 124 / 112 g n _|_ 112 Sn reactions in log-scale. From this 
figure, we observe that the ratio for each isotope Z exhibits a remarkable exponential 
behavior. For each isotope (Z), an exponential function form Cexp(a./V) is used to fit 
the calculated points and the parameters a are obtained for all isotopes. Analogous 
behavior is observed for each isotone (A), an exponential function form C exp(/3Z) is 
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used to fit the calculated points and the parameters (3 are obtained for all isotones. 
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Figure 4. Same as in Fig. [2] but for the isoscaling parameters a (upper panel) and 
\f3\ (lower panel) as a function of the excitation energy per nucleon. 

In Fig. HI we present the extracted slope parameters a (upper panel) and \(3\ (lower 
panel) of the exponential fits as a function of the excitation energy per nucleon (E*/A). 
Since the excitation energy changes as a function of time in the process of evaporation, 
the values of E* / A shown here are taken at the beginning of the evaporation stage and 
the mass of the prefragment is used to calculated E* / A. In the model the excitation 
energy is proportional to the abraded nucleons and can reflect the violence of the 
collision as the parameter (Zp ro j — Z) / Zp ro j shows [19]. In this figure, a and \(3\ 
show a decreasing trend with the increasing of E* / A. This behavior for projectile- 
like fragments is different with light particles. The isoscaling parameters of light 
fragments from multifragmentation is almost constant for different isotopes because the 
excitation energy or temperature is almost same for all light fragments in the process 
of multifragmentation. The decrease of the isoscaling parameters in our calculations 
may mainly be attributed to the evaporation effect of the prefragment with different 
excitation energy as in the disappearance of the isospin effect |T8], [19]. The values of a 
and |/3 1 are quite different from different reaction systems due to the different size and 
isospin of the projectiles. 

According to Eq. (2) and (3), a and (3 have a linear dependence on A[(-|) 2 ] 
or A[(-j) 2 ]. Since this parameter is dependent on the reaction system, we divide 
a (P) by A[(-|) 2 ] (A[(^) 2 ]) to remove the system isospin and size dependence and 
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Figure 5. Same as in Fig. [2] but for the reduced isoscaling parameters a/A[(^) 2 ] 
(upper panel) and /3/A[(^) 2 ] (lower panel) as a function of the excitation energy per 
nucleon. 



call them reduced isoscaling parameters. The results are given in Fig. [5j After the 
reduction, a/A[(^) 2 ] (/?/A[(^) 2 ]) of different reaction systems demonstrate almost same 
dependence with E* / A. Eq. (2) and (3) are deduced from the grand-canonical limit for 
multifragmentation of hot source. For projectile-like fragments, the same behavior is 
observed in the SAA model. In this sense, the reduced isoscaling parameters a/A[(^) 2 ] 
(/3/A[(^) 2 ]) may be used as a sensitive observable for measuring the excitation extent 
of projectile-like fragments during the collisions without system size dependence. 

From above discussions, we found that AApg^/ A Ap ro j and ct/A[(-|) 2 ] 
(/5/A[(^) 2 ]) decrease with (Zp ro j — Z) / Z^ ro j or the excitation energy per nucleon. But 
the later one decreases slower than the first one. It means that a/A[(-|) 2 ] (/3/A[(-^) 2 ]) 
is more sensitive to the isospin effect of the projectiles. Since a/A[(^) 2 ] (/3/A[(-j) 2 ]) 

is related to — j, — as in Eq. (2) and (3), it can be used as an observable to study the 
excitation and asymmetric nuclear equation of state in heavy ion collisions. 

If we use the Fermi-gas relationship between excitation energy per nucleon and 
temperature E* / A = ^T 2 to calculate T tentatively, with the inverse level density 
parameter a=8— 13 (in our calculation a=10 is used), the symmetry energy coefficient 
(Csym) could be extracted. Results extracted from a/A[(^) 2 ] and /3/A[(^) 2 ] are shown 
in Fig. [61 For E* / A around 1 MeV, the symmetry energy coefficients from a and (5 are 
around 15 MeV. These values are a little lower than the standard value Csym=25 MeV 
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E*/A (MeV) 

Figure 6. Same as in Fig. [2] but for the extracted symmetry energy coefficient from a 
(Cgymi upper panel) and j3 (Cgy m , lower panel) as a function of the excitation energy 
per nuclcon. 



in liquid drop model [23], but it seems consistent with the extracted results from the 
experimental data by A. Le Fevre et al. [28]. The obtained Csym is not a constant 
and decreases with the increase of E* / A. Similar dependence was also observed in 
the experimental studies but their Csym values is a little bit larger than ours [28, 29J. 
It should be pointed out that the experimental data are taken at different incident 
energies (300 A MeV and 600 A MeV in Ref. [28], around 25 A MeV in Ref. [29]). 
Our calculations are performed at 60 A MeV, but we have found that there is almost 
no incident energy dependence for the isoscaling parameters in our model. In our 
results, the symmetry energy coefficient decreases quickly when E* / A is less than 1 
MeV. This may stem from the increase of the inverse level density parameter at low 
excitation energy [30]. Of course, in this low E* j A range there is very few theoretical 
and experimental data up to now and more researches are necessary. Experimentally it 
may be difficult to extract E* / A. As we have mentioned previously that the parameter 
(Zp ro j — Z)/Zp ro j can reflect the violence of nuclear collision and is approximately 
proportional to E*/A in not very central collisions. Actually quite similar dependence 
as in Fig. 0H6]is seem if (Zp ro j — Z)/Z^ TO j is used instead of E*/A. Thus we can also 
study the dependence of the isoscaling parameter and symmetry energy coefficient with 
(Zp ro j — Z)/Zp ro j experimentally when there is no E*/A data. 

However, some cautions should be reminded. In the SAA model, the symmetry 
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energy term is not taken into account explicitly. But a similar analysis of isoscaling as 
in the statistical multifragmentation model could be done since there exists different 
isotopic and isotonic distributions between two systems. Of course, the effect of 
symmetry energy term should have been reflected implicitly in the assumption of the 
abrasion and also evaporation stages in the SAA model. Thus a clear isoscaling behavior 
is observed for projectile-like fragments in the this work and the symmetry energy 
coefficient extracted based on Eq. (2) and (3) is consistent with the experimental 
data [281 [29]. The present calculation could provide some useful information for further 
experimental and theoretical investigations on the isoscaling of projectile-like fragments. 



4. Summary 

In summary, systematic behaviors of the isospin effect and isoscaling of projectile-like 
fragments from 40 / 36 Ar, 48 / 40 Ca, 64 / 58 Ni, 86 / 78 Kr, 124 / 112 Sn and 129 / 136 Xe on 112 Sn at 
60A MeV have been studied by a modified statistical abrasion-ablation model. The 
normalized peak differences AAp ea ^ / AAp ro j for different reaction systems show similar 
dependence with the parameter (Zp ro j — Z) j Z-p TO y The isoscaling parameters a and j3 
are extracted for the produced isotopes and isotones, and they show different values for 
different systems. However, the reduced isoscaling parameters a/A[(^) 2 ] and /5/A[(^) 2 ] 
show almost same dependence with E* / A for different systems. Assuming a Fermi-gas 
behavior, the symmetry energy coefficients are tentatively extracted from a and f3 and 
it seems that the results are consistent with the experimental data. 
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